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ABSTRACT 

We study spectroscopically determined iron abundances of 642 solar-type stars to 
search for the signature of accreted iron-rich material. We find that the metallicity 
[Fe/H] of a subset of 466 main sequence stars, when plotted as a function of stellar 
mass, mimics the pattern seen in lithium abundances in open clusters. Using Monte 
Carlo models we find that, on average, these stars have accreted ~ OAM^ of iron while 
on the main sequence. A much smaller sample of 19 stars in the Hertzsprung gap, 
which are slightly evolved and whose convection zones are significantly more massive, 
have lower average [Fe/H], and their metallicity shows no clear variation with stellar 
mass. These findings suggest that terrestrial-type material is common around solar type 
stars. 

Subject headings: planetary systems — stars: abundances — stars: chemically peculiar 



1. INTRODUCTION 

The discovery that at least 6 — 8% of solar type stars harbor Jupiter-mass or larger bodies, of- 
ten in small, eccentric orbits (Mayor & Queloz 1995; Marcy & Butler 1996; Butler & Marcy 1996; 
Marcy et al. 2000) shows that planets are not exceptionally rare. The transiting planet orbiting 
HD 209458 has a mass of 0.69 Jupiter masses, and a radius about 1.4 times that of Jupiter. Clearly 
it is a gas giant like Jupiter or Saturn; the minimum masses of the other known objects suggest 
that they are also gas giants. The Doppler technique used for most of the discoveries cannot find 
terrestrial mass objects in AU scale orbits. What fraction of solar type stars have terrestrial- type 
planets? 
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How could an observer with our current technology located on a star in the solar neighborhood 

decide if there were any terrestrial-type bodies orbiting the sun? More generally, how could such 
an observer estimate the fraction of solar type stars having companions of terrestrial (as opposed 
to gas giant) compositions? In this paper we explore one possible way of addressing the latter 
question, stellar-mass dependent photospheric metallicity trends. 

In section 2 we argue that a few Earth masses of rocky /icy material have accreted onto the 
Sun over its lifetime. In section 3 we discuss one possible observational tracer of similar accretion 
occurring on other solar type stars. In section 4 we examine a sample of stars with spectroscopically 
measured photospheric iron abundances 

(1) 

where fpe is the mass abundance of iron in the photosphere of the star, and /fb,© ~ 1.3 x 10~^ 
is the the mass abundance of iron in the photosphere of the sun. We compare the observations to 
Monte Carlo models of stellar pollution in section 5. We discuss the implications of our results in 
section 6 and present our conclusion in the final section. 

2. POLLUTION BY SMALL BODIES IN OUR SOLAR SYSTEM 

The surface density of iron in our solar system follows a rough power law with distance 
from the sun (Wcidcnschilling 1977). We have compiled recent estimates for the iron content of 
the planets in table 1. From a least squares fit including Venus, Earth, and the four gas giants we 
find Si?e ~ 4{r /I AU)"^'"^ gcm'~^ . Weidenschilling noted a clear deficit of material in the region 
interior to Venus, and between Earth and Jupiter, relative to this power law. Gas drag acting on 
planetesimals near the sun could reduce the surface density of such bodies near the present orbit 
of Mercury by dropping them onto the sun. Because the solar convection zone was very deep for 
~ 20 million years (see figure 1), while the gas disk that would have produced the drag is believed 
to have survived only a few million, this would not have altered the apparent metallicity of the sun. 

Materiel in the asteroid belt clearly survived such gas drag. However, Kirkwood noticed 
over one hundred years ago that the distribution of asteroids showed distinct gaps at the loca- 
tion of orbital resonances with Jupiter (Kirkwood 1867). This suggests that material has been 
removed from the gaps under the influence of Jupiter. Over the last twenty years, the dynamics of 
this removal have been worked out in considerable detail (Wisdom 1983; Holman &; Murray 1996; 
Gladman et al. 1997). The gaps are the result of resonant, chaotic perturbations of the asteroid 
orbits by Jupiter. A second feature, seen in plots of orbital eccentricity versus semimajor axis, is 
the result of a secular resonance, the i'q resonance. In either type of resonance bodies with semi- 
major axis larger than about half that of Jupiter's tend to be removed from the solar system; their 
eccentricities grow with time, increasing their apoapses until they reach the orbit of Jupiter. Close 
encounters with that planet then quickly eject the asteroids. 



[Fe/H] = log 



fPe 
jFe,e. 
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Resonant asteroids in smaller orbits tend not to be ejected; their apoapse does not reach 
the orbit of Jupiter. Instead, their periapscs decrease until they hit the sun. Several related 
mechanisms can combine to produce this decrease in periapse, but the two main actors are mean 
motion resonances and the resonance already mentioned. In a mean motion resonance such as 
that responsible for the 3/1 Kirkwood gap, the eccentricity of the resonant asteroid undergoes a 
random walk. Since the eccentricity cannot decrease below zero, there is a tendency for bodies 
with small e to drift to larger e. When the eccentricity is large enough, the asteroid strikes the 
sun (Gladman et al. 1997). A small fraction will suffer close encounters with the Earth or other 
terrestrial bodies; this may lead to Jupiter crossing orbits, which usually leads to ejection of the 
asteroid from the solar system. Only a very small fraction will collide with either a terrestrial planet 
or Jupiter. In the secular resonance the precession rate equals the sixth secular frequency of 
the solar system. This frequency appears with substantial amplitude in the orbital motion of both 
Jupiter and Saturn; it is roughly the average precession rate of the latter planet. The resonance 
rapidly removes angular momentum from the orbit of the asteroid, dropping the periapse below 
the surface of the sun. 

These resonances are currently active, and arc believed to be responsible for the bulk of the 
meteorites striking the Earth. While there is no direct observational confirmation, it is virtually 
certain that a much larger flux of asteroidal material from these resonances is currently hitting the 
sun. 

It is also likely that the flux of asteroidal material was much larger in the past. The locations 
of both mean motion and secular resonances depends on the positions of Jupiter and Saturn, and 
on any massive gas disk. It has been suggested that as the protoplanetary gas disk dissipated 
the 1/6 resonance swept through the belt, depleting much of the mass. More recently it has been 
suggested that one or both of Jupiter and Saturn migrated to their current positions, causing 
a similar sweeping of mean motion resonances. This would preferentially deplete the outer belt 
(Holman &; Murray 1996; Liou & Malhotra 1997), but would also drop material from the inner 
belt onto the sun. Finally, we note that the asteroid belt is hot, in the sense that the distribution of 
both e and i is broad. A number of authors have suggested that a planet size body swept through 
the region early in the history of the solar system, heating the belt as well as dropping material on 
the star (Wetherill 1992; Petit et al. 1991). 

The surfaces of the terrestrial planets and the Moon show evidence for an extended period 
of bombardment, known as the late heavy bombardment, lasting up to a billion years after their 
surfaces cooled. This is consistent with some of the longer lived depletion scenarios described above. 

From table 1, the amount of material inferred to have populated the original asteroid belt is 
of order five to ten Earth masses. Roughly 2-5 Earth masses would have been between the present 
position of Mars and 2.5 AU, half the present semimajor axis of Jupiter. A substantial fraction of 
this material, of order half, would have ended up in the sun. We will take ~ 2M^ as a representative 
number. Meteoritic material is roughly 20% iron by weight (Grevesse & Anders 1989), so we take 
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~ O.5M0 as the amount of iron that was dropped on the sun after the convection zone thinned. 

3. POSSIBLE OBSERVATIONAL EVIDENCE FOR POLLUTION 

There is another possible piece of evidence besides the late heavy bombardment for such a 
late depletion of the asteroid belt. As we noted above, a substantial fraction of the "missing" 
material originally in the asteroid belt strikes the sun. This material will be mixed throughout the 
convection zone of the star. Since a substantial bombardment of the terrestrial planets lasted for 
much longer than 100 million years, this material would have landed in a convection zone of mass 
~ 3 X 1O~^M0. Hence there is a possibility that the signature of solar bombardment was left in 
the form of an enhanced metallicity in the solar envelope. 

The mass of the convection zone was ~ 3% of the solar mass between 10^ and 10^ years, while 
the present mass fraction of iron observed in the photosphere is 1.3 X 10~^. This yields a total 
iron content of ~ 12.75M0 between 10^ and 10^ years. Adding half an Earth mass of iron would 
increase the observed [Fe/H] by only « 0.017 dex. 

Precision measurements of the bulk metallicity of a single solar type star are not presently 
available. Measurements of the sound speed in the interior of the sun, using the five minute 
oscillations, could in principle provide such a measurement. The only work along these lines that 
we are aware of is the paper by Henney & Ulrich (1998). They show that the addition of ~ SM® 
of cometary material results in a frequency shift that, because of uncertainties in solar models, is 
too small to be detected. 

However, it is possible that observations of large numbers of stars could reveal the presence 
of pollution. For example, stars having a fixed metallicity but different masses will mix material 
dropped onto their surfaces to different depths. Knowledge of the depth of this surface mixing 
layer, which we summarize in the next subsection, allows us to look for variations of metallicity 
with stellar mass produced by accreted material. 

3.1. The Mass Of The Surface Mixing Region 

Material dropped onto the surface of a star will not in general be confined to the photosphere. 
The surface layers of stars less than about IAMq are convectively unstable, and any material 
dropped on the star is expected be mixed throughout the convection zone. Observations of lithium 
(Li) abundance provide support for this theoretical picture. But these same observations suggest 
that there is another form of mixing in stars more massive than ~ 1.2Mq, possibly related to 
meridional circulation. 

The evidence for the two types of mixing comes in the form of variations in the surface lithium 
abundance with stellar mass. Lithium is destroyed when it is exposed to temperatures above about 
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3.1 X lO^K; the temperature at the photosphere is much lower than this, so lithium depletion 

indicates that photospheric material is mixed down to regions where the temperature is higher. 
Lithium could also be removed from the photosphere by simple settling (so that the lithium is 
stored below the photosphere, as opposed to being destroyed) or by stellar mass loss. However, the 
observations, which we review below, point to destruction by mixing. 

Observations show that for stars with mass below ~ 1.2Mq the photospheric abundance of 
lithium decreases rapidly with decreasing stellar mass, e.g., (Boesgaard 1991). Some depletion 
occurs before the main sequence, but there is strong evidence that depletion also occurs on the 
main sequence (Jones et al. 1999). The main sequence depletion is believed to be related to the 
increasing depth of the bottom of the surface convection zone with decreasing stellar mass in these 
stars, supplemented by some convective overshoot and/or settling in stars near the upper end of 
this mass range. 

Stars above about 1.2Mq have very thin, or even lack, convection zones. In standard stellar 
models, with no settling, mixing or mass loss, the surface Li and Be abundances are constant 
over the main sequence lifetime of such high mass stars. Both Li and Be are destroyed below 
modest depths in the interiors of such stars. For stars in the mass range 1.2-1.6Mq and having 
solar metallicity, our stellar models indicate that the nominally undepleted region has a mass of 
> 3 X 10~^Mq for Li and > 6 x 10~'^Mq for Be. (We assume that the lithium is destroyed at 
temperatures above 3.1 x lO^K). At depths below those corresponding to these masses, Li and Be 
should be destroyed. 

Real stars with masses larger than 1.2Mq don't behave like the models. Hyadcs stars show 
a distinct "lithium dip" centered around IAMq (Boesgaard &; Tripicco 1986). Between ~ 1.2 and 
IAMq the lithium abundance is seen to decrease with increasing stellar mass. Prom 1.4 to 1.5Mq 
there is a sharp increase in lithium abundance. Balachandran (1995) finds the highest abundances 
for stars with masses above 1.5Mq (the "blue side" of the dip) and for stars of 1.2M0 (the 
"red side" of the dip); these stars have abundances about equal to those found in meteorites 
(Grevesse & Anders 1989), abundances which are generally believed to be primordial. 

Observations of young stellar clusters such as a Per and the Pleiades find little evidence 
of lithium depletion in more massive stars (Soderblom et al. 1993). Indeed, there is a strong 
correlation between cluster age and Li abundance in the "lithium dip" between 1.2 and 1.5Mq, 
indicating that the depletion takes place on the main sequence on a time scale of 100 Myrs or more 
(Boesgaard 1991). 

Observations of clusters having different metallicities [Fe/H] show that the location of the dip is 
slightly metallicity dependent (Balachandran 1995). Lower metallicity clusters have dips centered 
at lower masses. However, the high mass side of the dip depends only weakly on metallicity. 

As we noted above, the lithium depletion appears to be due to mixing rather than gravitational 
settling or stellar mass loss. Two observations lead to this conclusion. First, surface depletions of 
lithium and beryllium are correlated, with the surface lithium diminishing more rapidly than surface 
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beryllium (Deliyannis et al.). Because the Be rich zone has twice the mass of the Li rich zone, mass 
loss would deplete all the Li before depleting the surface Be. Similarly, settling predicts that surface 
Li and Be would be depleted at roughly the same rate. Second, evolved stars with 1.2 < M < 1.5 
are also deficient in lithium (Gilroy 1989; Balachandran 1995); settling of lithium would leave a 
non-depleted layer of lithium below the surface (in the absence of mixing), which would be dredged 
up to the surface when the convection zone deepens. The fact that the predicted increase in evolved 
stars is not seen indicates that the Li missing from main sequence stellar photospheres is destroyed 
rather than sequestered below the photosphere. 

Above I.SMq there appears to be little mixing of the upper layers of the star with deeper 
layers, since there is little Li depletion. However, the lithium in these stars does appear to be 
destroyed in deeper layers; evolved (giant) 1.6 solar mass stars are severely depleted (Gilroy 1989). 
This depletion occurs inside the progenitors, below the photosphere (Vauclair 1991), so that when 
the surface convection zone forms and then thickens, the surface lithium is mixed over a larger 
volume of depleted material, reducing the photospheric lithium abundance. 

The surface mixing seen in the dip appears to be rotationally induced (Charbonnel et al. 1992; 
Balachandran 1995), although this is far from certain. 

As already noted, stars in a Per and the Pleiades clusters, with estimated ages of 50 and 
70 million years, show little depletion in Lithium abundance for stars between 1.2 and 1.5 solar 
masses, while stars in the Hyades and older clusters show depletions that increase with cluster age. 
We interpret this to mean that the upper layers of the stars are mixed downward on time scales of 
^ 100 million years. 

The observational data summarized here leads us to propose the following simple model for 
the mass M^ix of the surface mixing region in solar mass stars (see Figure 2). We first describe 
the lithium dip. The location, in mass, of the dip is described by Mred, Mcentral, and Muue, giving 
the low, central, and high mass points of the dip. These three masses depend on the stellar [Fe/H]. 

Above Alblue we assume M„iix = M^^; we use M^/i = 3 x IO^^Mq, the value obtained for the 
meridonal circulation model of Charbonnel et al (1992). This number is currently unconstrained 
by observation on the low side, but must be less than !=a 5 x lO^^M© to prevent Li depletion. 

For Mcentral ^ < we use a linear interpolation between M^/j and M^c ~ 3 x 10~^, 

where the latter value (for a star at the center of the dip) allows for complete Li depletion but only 
partial Be depletion. We also use a linear interpolation between M^nc and M^i for M^g^ < M* < 
Mcentral- Stars with M = Mred appear to be only slightly depleted, so we take Mj^i ~ 5 x lO^^M©. 

Extrapolating from data in (Balachandran 1995), we take 



Mred/MQ 
Mcentral I Mq 
Muuc/Mq 



1.51 + 0.054[Fe/H] 
1.42 + 0.25[Fe/H] 
1.31 + 0.31[Fe/H]. 



(2) 
(3) 
(4) 
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Note that the high mass or blue side of the dip is essentially independent of [Fe/H]. For stars with 
mass below Mredi we take the mass of the convection zone (Mcvz) plus the mass of the mixing zone 
Mmi of a Mred Star on the low mass or red edge of the dip. 

The mass of the convection zone at an age of 10® years (the mass at 10^ is similar, see Figure 
1) is found from the following empirical fit, derived from our stellar models: 

Mcv, = a ([Fe/H]) x (M_3 - M,)^•6^ (5) 

where a ([Fe/H]) = -0.2[Fe/H] + 0.49. The quantity M_3 is the (metalhcity dependent) stellar 
mass at which the convection zone mass equals IO^^Mq, 

M_3 = a_3[Fe/H]+6_3, (6) 

where a_3 = 0.32 and 6_3 = 1.45. Equation (5) is very accurate near M_3 and good to about 20% 
at O.SMq. 

Our model for the depth of the mixing region in a solar metallicity star is shown in Figure (2). 
Note that stars of different metallicity will have "lithium dips" (peaks in this figure) at slightly 
different masses than those shown in the figure. 

Armed with this model, we can predict the average change in metallicity of a sample of stars 
subjected to iron accretion. Essentially we will look for behavior similar to that seen in the abun- 
dance of lithium. Polluted low mass stars of a fixed age will show an increase in metallicity with 
increasing mass (below the low mass end of the Li dip, around 1.0 solar masses). Intermediate 
mass polluted stars will show a decrease in [Fe/H] with increasing mass, up to M^ip, followed by a 
sharp increase at Muue- 

We want to stress that, while we are using the lithium data as a proxy to estimate the depth 
of the surface mixing layer of main sequence solar type stars, and while we anticipate that the 
surface abundances of lithium and iron will follow similar patterns, the reason for the abundance 
variations differ. The lithium abundances vary because lithimum is destroyed when it is mixed 
below the surface of the star. Stars in the mass range we consider neither destroy nor create iron. 
Rather, we assume that some iron is added to the surface of the star after it forms; it is then diluted 
by the mixing process down to the same depth that the lithium is mixed to. In other words the 
variations in lithium are produced by the destruction of lithium inside the star, while the variations 
in iron are produced by the addition of iron from outside the star. 



4. OBSERVATIONAL CONSTRAINTS ON POLLUTION 

In order to look for mass-dependent variations in [Fe/H], we must be able to estimate the 
masses of stars, and we must have reliable measurements of the metallicity. To find the mass of a 
star we need to know its luminosity, or alternately its absolute V magnitude My, its colors (we use 
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B—V), and its composition [Fe/H]. The HIPPARCOS catalogue gives parallaxes accurate to roughly 

one milliarcsecond for a selection of stars in the solar neighborhood, allowing us to find absolute 
magnitudes when we are given V. We use SIMBAD values for B and V. We take metallicities 
from the (Cayrel de Strobel et al. 1997) catalogue. This lists spectroscopically determined values 
for [Fe/H] taken from the literature. 

We take all the HD stars in the Cayrel de Strobel et al. catalogue for which there are HIPPAR- 
COS parallaxes larger than 10 milliarcseconds (corresponding to distances less than lOOpc). We 
also require that the error in the parallax be less than 10 percent. Since we are interested in main 
sequence stars that have masses less than about 2Mq, we eliminate all stars with My < 1.0. We 
then fit the remaining stars to stellar tracks taken from a grid of models having —0.5 < [Fe/H] < 0.5 
in steps of 0.05 dex, and masses in the range 0.60 < M < 1.75Mq in steps of O.O5M0. The models 
are from the Yale stellar evolution code. 

A star is defined to be evolved if it has a convection zone mass more than ten times larger 
than the mass of the surface mixing layer at an age of 10*^ years, while the absolute V magnitude 
has not decreased by more than pa 0.5. We refer to such stars as "Hertzsprung gap" stars. This 
occurs roughly when the star leaves the main sequence. It is a useful definition because it allows us 
another check on the pollution scenario; if pollution is occurring, the class of stars with convection 
zone masses ten times larger than their initial main sequence convection zone masses should show 
substantially lower [Fe/H] values than their parent populations. 

If the absolute V magnitude has decreased by more than pa 0.5, the star is considered to be 
a giant. Because their surface gravities are so much lower than main sequence stars of the same 
mass, we hesitate to compare the metallicity trends of giants and dwarfs. 

In the course of our analysis a problem with our stellar models became apparent. Stars with 
B — V > I were systematically found to be older than 20Gyrs. We believe that this is due to a 
failure of our atmospheric models at these low effective temperature. For this reason we do not 
include any stars redder than this in our analysis; this corresponds to stars with masses less than 
~ 0.75 — O.8OM0. At the other end of the mass spectrum, we felt confident in extrapolating up to 
masses of I.8OM0, O.O5M0 beyond our grid. 

Of the 650 stars in our sample, 466 were unevolved and had masses in the range 0.8 < M < 1.8. 
Fifteen stars had convection zones 3-10 times larger than their main sequence value (at lO^yrs); we 
treated these as unevolved. Nineteen stars had convection zones more than ten times as massive 
as the main sequence values; these are our Hertsprung gap stars. 

In Figure (3) we plot [Fe/H] against stellar age, while Figure (4) shows [Fe/H] plotted against 
stellar mass. The metallicity decreases with increasing age, and increases with increasing mass. We 
also note that the Hertzsprung gap stars tend to have lower metallicity than unevolved stars of the 
same age and mass. Figure (5) shows the histogram of [Fc/H] for our sample. It is well-fit by a 
gaussian with mean ([Fe/H]) = —0.095. 
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The trend with stellar age is expected on the basis of chemical evolution; younger stars are 
constructed using gas that has been contaminated by material processed in the interiors of massive 
stars. The more recent the stellar birth, the higher the fraction of processed material. The trend 
is seen even more clearly when we bin the data by age, as seen in Figure (6). The slope of a least 
squares fit to the [Fe/H] versus log(age) curve is —0.21. 

The trend of increased metallicity with increased stellar mass might arise from the age- 
metallicity trend, combined with the short lifetimes of massive stars. Since the stars we consider 
are either unevolved or only slightly evolved, they must have ages less than or roughly equal to 
their main sequence lifetimes. Massive stars have short lifetimes, and hence low ages (see Figure 7). 
They will therefore be more metal rich, on average, than less massive stars. 

In the next section we will show that, while part of the variation of metallicity with stellar 
mass is due to the age-metallicity relation, it is not the whole story. 



The data suggest a possible trend of surface [Fe/H] with stellar mass. In order to evaluate the 
significance of this trend, we have undertaken Monte Carlo experiments to calculate the predicted 
metallicity of a population of polluted stars. 

We start by producing a population of unpolluted stars having the same characteristics as the 
Cayrel dc Strobcl-dcrivcd sample described above. This involves generating a sample of stars with 
mass roughly uniformly distributed between 0.8 and I.SMq, with mctallicities gaussian distributed 
about an age dependent mean. Prom our sample of stars and using a simple two parameter fit we 
find that the mean [Fe/H] is given by 



where a = — 0.21zb0.02 and b = 0.03. The width of the [Fe/H] distribution is slightly age dependent, 
but we assume an age-independent width a = 0.18. 

We constrain the stellar age to be less than the calculated lifetime of the star. As noted above 

we define the lifetime of a star to be the age at which the mass of the surface convection zone 
exceeds ten times the mass of the mixed surface layer at an age of 100 Myrs. Our stellar models 
have lifetimes (as defined above) which vary as the 3.65 power of the stellar mass. 



5. 



MONTE CARLO MODELS OF POLLUTION 



{[Fe/H]) 



a X log{Age) + /?, 



(7) 





where 



Mi([Fe/H]) = 11.5 + 4.5[Fe/H]Gyr 
is the lifetime of a one solar mass star of metallicity [Fe/H]. 



(9) 
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using this prescription we generate 466 stars, yielding the mass-age relationship depicted by 
the open circles in Figure 7; the agreement with the observed age-mass relationship is very good. 



We can use this Monte Carlo model to see if the [Fe/H]-mass correlation seen in our sample 
could be due simply to the known age-mass and age-[Fe/H] correlation. To repeat, the idea is as 
follows: massive stars are necessarily young. Young stars are metal rich, since they form out of 
gas that has been processed through many generations of high-mass (supernova producing) stars. 
Thus more massive stars should have higher metallicity than less massive stars. 

This argument is partially borne out by the Monte Carlo model. We have adjusted the slope 
of the [Fe/H] vs. logarithm of stellar age relation to obtain agreement with that seen in our sample 
(—0.21 dex per log(Gyr)), then examined the metallicity as a function of stellar mass. The stellar 
metallicities do show a dependence on stellar mass. However, the slope of the [Fe/H] vs mass 
relation for the observed stars is 0.26 it 0.03 dex per solar mass, while that of the unpolluted Monte 
Carlo model is 0.18 ± 0.03, different by 3 standard deviations. 

A second indication that there is a real mass dependence in the data comes from allowing 
for a third parameter in a least squares fit. Heretofore we have used a two parameter fit, an age- 
metallicity slope and an intercept. If we allow for a linear variation of [Fe/H] with mass, the slope 
of the fit is a third parameter. A fit of the form 



does give a significantly improved chi squared. 

To see if this three parameter fit can really distinguish between a pure age dependence and 
a mixed age and mass dependence, wc modified our Monte Carlo model to allow for a linear 
dependence of average metallicity on stellar mass. We then ran the resulting stars through our 
least squares routine. If we set the amount of added iron to zero but force an age dependence, 
the least squares fit finds the correct age-metallicity slope. It also correctly finds that there is 
no mass dependence, despite the apparent dependence in a metallicity versus mass plot. If we 
introduce a mass dependence into the model, the least squares fit correctly finds both the age 
and mass dependences. Finally, we produced samples of stars with a mass dependence but no age 
dependence. A plot of metallicity versus age clearly shows a trend; this is expected since we force 
the more massive stars, which are younger, to have higher metallicities. Despite this apparent 
age-metallicity trend, the linear least squares fit finds the correct mass-metallicity slope, and a zero 
slope for the age-metallicity relation. 

These results suggest that the real data exhibits a mass trend independent of stellar age, though 
by itself the argument is not very convincing. More compelling evidence is obtained by comparing 



5.1. The Unpolluted [Fe/H]-Mass Correlation 




(10) 
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to more sophisticated models. The inadequacy of a simple age dependent mctallicity is most clearly 
seen by examining the binned metallicity versus mass distribution, Figure (8). Another indication 
is the poor /dof of the fit; the best unpolluted model, shown as the dotted line in the figure, 
gives = 3.2 for the binned distribution. We proceed to examine more realistic polluted models, 
which give substantially better fits to the data. 

5.2. Realistic Polluted Models 

As noted above, pollution introduces a third free parameter, the amount of iron added to the 
star. In the last subsection we assumed that [Fe/H] increases linearly with stellar mass. This is 
not what we expect from pollution, since the mass of the surface mixing layer is not linear, or even 
monotonic with stellar mass. In this section we assume a gaussian distribution of accreted mass, 
with a variance equal to half the mean. Stars allotted a negative amount of accreted mass by this 
process are assumed to have no added material. We assume the accreted mass is mixed over a 
surface layer of mass M^^^, as described in section 3.1. Our model for the mass of the surface 
mixing layer has many parameters, but we hold them fixed at the values given in that section. The 
single adjustable parameter is then the mean mass of accreted iron. We adjusted this mean to 
obtain a minimal in the binned data, both metallicity-mass and mctallicity-agc. The resulting 
fit to the mctallicity- mass data is shown as the solid line in Figure (8). We were able to find fits 
with reduced equal to one. 

The reduction in the of the fit afforded by the more sophisticated pollution model is dra- 
matic. The detailed variation of [Fe/H] with stellar mass, including the steep increase starting at 
O.SMq, the dip around IAMq, and the steep increase to l.SM© seen in the binned data is rea- 
sonably similar to that seen in the model. We feel that the low x^ ^^'^ the detailed variations of 
metallicity with mass offer strong support both for 1) a mass dependent mixing, as suggested by 
the Li dip; and 2) for accretion of ^ 0.5 ± 0.2M^ of iron in most stars in the solar neighborhood. 

5.3. Hertzsprung-gap Stars 

The metallicity-mass trend, combined with the jump in average [Fe/H] around 1.5 solar masses 
strongly suggests that stars in our sample have accreted iron rich material. However, the properties 
of our sample are not well constrained, since the selection criteria of the various samples collated by 
Cayrel de Strobel are difficult to ascertain. Perhaps the apparent variations in average metallicity 
are due to some unknown selection effect. For example, one of the subsamples is that of Favata, 
Micela, & Sciortino (1997). A plot of [Fe/H] versus mass for this sample shows an increase toward 
large masses, starting at about 1.15M0. However, this sample was selected to have B — V > 0.5; 
this eliminates low metallicity high mass stars, which are bluer than this limit. This is the origin 
of the increase in [Fe/H] at 1.15M0 in this subsample. 
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We do not see an increase in [Fe/H] at I.ISM© in the large sample. However, perhaps the sharp 
jump we see at 1.5Mq is due to a similar selection effect. The best way to eliminate this possibility 
would be to perform an unbiased survey of a large number of stars in the solar neighborhood. We 
encourage observers to undertake such a survey. 

We do have available to us a small control sample. Stars that are slightly evolved, in the 
sense that their convection zones are 10 or more times larger than those on the main sequence, but 
which are not yet giants (having My within 0.5 of the main sequence value) provide a check on 
our interpretation. These stars sit in the Hertzsprung gap. Unlike giants, their surface gravities 
and fluxes are very similar to stars on the main sequence, so their metallicities can be compared 
directly to those of main sequence stars. 

If pollution is responsible for the [Fe/H] variations over and above those due to age effects, 
than Hertzsprung-gap stars should have lower average metallicities, and they should show no trends 
with stellar mass. 

The open triangles in Figure (4) represent the metallicities of stars in the gap. It is apparent 
that these stars have, on average, slightly lower metallicity than unevolved stars of similar mass. 
This is confirmed by an examination of the binned data. Figure (8), which also shows that there 
is no evidence for a change in average [Fe/H] in the gap stars at the blue edge of the lithium gap 
around 1.5 solar masses. This provides support for the notion that the variations in metallicity 
with mass seen in the unevolved stars is due to pollution. 

The dashed line in Fig. (8) shows the predicted [Fe/H] for the same model as the solid line, 
but where the mass of the accreted material has been reduced by a factor of ten. This is equivalent 

to increasing the mass of the surface mixing layer by the same factor. The predicted metallicity 
for this population is lower than that predicted for the polluted sample, and shows only a hint of 
a "lithium" dip. It gives an acceptable fit to the rather sparse observational data (reduced less 
than one). 

In plotting both the data and the fit we have not adjusted [Fe/H] according to the star's age. 
In Figure (9) we replot the data adjusting for the age- metallicity trend using the measured slope 
of the age-metallicity relation in Figure (6). This involves increasing the [Fe/H] for old stars. This 
figure should be compared to plots of lithium abundance in clusters; we see evidence of a "lithium" 
dip in the iron data similar to that seen in lithium data. The dip is less distinct in the iron data, a 
fact we attribute to the range of metallicities in our sample (the metallicities in cluster stars show 
much smaller dispersions). The variation in metallicity produces a variation in the mass of the 
dip, so combining stars with different metallicity will tend to smear out the dip. This smearing is 
clearly seen in our Monte Carlo models (the solid line in Fig. 8). 

We have also examine the metallicities of giants in our sample. They show a slight increase 
in [Fe/H] with mass, but this increase is consistent with the age-[Fe/H] correlation. The average 
[Fe/H] is actually slightly larger than either the unevolved stars or the Hertzsprung gap stars, 
although consistent within the rather large errors. However, given the very different environment 
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for line formation in these stars, we do not feel it is appropriate to make a direct comparison 
between the two populations. 

6. DISCUSSION 

We made use of three free parameters in fitting the metallicity-agc and metallicity-mass data; 
one for the slope of the metallicity-age relation, one (the mean of the added mass) for the correlation 
in the metallicity-mass relation, and an overall normalization. We want to stress that, in fitting 
the Hertzsprung gap star data, (both age and metallicity) we did not adjust any free parameters. 

We also had a number of parameters in the model which were fixed by other observations; 
most were used to describe the location and depth of the lithium dip. The most critical for our 
purposes is the mass M^nh = 3 x 10~^Mq of the surface mixing layer in stars blueward of the dip. 
Varying this mass will directly affect our estimates of the amount of iron accreted onto the star. 
For this reason our estimate of OAMq of accreted iron is uncertain. However, the need to accrete 
some material is not uncertain; reducing M^h will reduce the inferred mass of accreted material, 
but models with no pollution will still fail to fit the data. 

The linear model for the mass of the surface mixing layer we employ is unlikely to be correct 
in detail; the "v" shaped bottom is particularly questionable. Future work should employ a more 
realistic model for Mmix, but for this first foray we did not feel that more realistic models were 
necessary. 

If the trends we see are due to pollution, one might surmise that stars with larger [Fe/H] 
will form with more massive planetesimal disks. This might result in larger amounts of iron rich 
material falling on stars with larger intrinsic metallicities. This could be modeled by assuming a 
correlation between [Fe/H] and the amount of added material; once again we leave such experiments 
for later work. 

6.1. Gas Migration and Ingestion of Gas Giant Planets 

The gas giant planets in our solar system appear to have rock/ice cores with masses of order 
IOMq, with the possible exception of Jupiter (Guillot 1999a), which could lack such a core. These 
cores are believed to form at > 5 AU from the sun since the surface density of solid matter in 
the protoplanetary disk is believed to jump up beyond that distance due to the presence of ice. 
The discovery of Jupiter-mass planets in small (less than 0.05 AU) orbits around solar type stars 
motivated a number of groups to suggest that gas giant planets migrate. If this migration is overly 
efficient, Jupiter-mass objects may be accreted onto the star. In the leading theory of planetary 
migration, that of tidal interactions between massive planets and the gas disk out of which they 
formed, the migration time is given by the viscous evolution time of the disk, prompting many 
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authors to predict such accretion (Lin 1997; Laughlin & Adams 1997). Others have suggested 
accretion of iron rich material originally between the planet and the star (Gonzalez 1997). 

There are two distinct scenarios. In the first, iron rich material is pushed onto the star by the 
gaseous disk. This material necessarily falls on the star before the gas disk disappears, when the star 
is between one and ten million years of age. In the second, the enriched material arrives via some 
other mechanism, possibly long after the gas disk disappears. The primary difference between the 
two, from the point of view of pollution, is the depth of the stellar convection zone; stars younger 
than ten million years and less massive than ~ 1.3 — IAMq have very massive convection zones 
and consequently are difficult to pollute. 

However, as pointed out by Laughlin and Adams (1997), stars more massive than 1.3 — IAMq 
will show evidence for pollution if they accrete iron rich material at ages ~ 10 Myrs. Figure (10) 
shows a Monte Carlo model where lAM^ of iron is added to stars before they reach ten million 
years of age. Wc have adjusted the age-metallicity slope and the fraction of stars accreting material 
to achieve the best ^t; we can find acceptable fits (with ~ 1) to the unevolved stellar data 
as illustrated by the solid line in the figure. Roughly 1/4 of these massive stars must accrete 
Jupiter-analogues as they approach the main sequence. 

This early accretion model does not do so well when compared to the evolved stars. The 
model predicts that the metallicities of the massive stars will lie below those of unevolved stars, as 
observed. However the model predicts that low mass stars will not show any reduction in ([Fe/H]), 
nor is the drop predicted for the massive stars as large as that seen in our data. Using our best fit 
model (the solid line in Figure 10) and then increasing the mass of the surface mixing layer by a 
factor of ten gives the dashed line in the figure; for this model the reduced = 2.5; the predicted 
metallicities are systematically higher than those observed. 

There is also a theoretical difficulty with the early accretion model. The model is motivated 
by the possibility that gas in the protoplanetary accretion disk can push Jupiter- mass objects into 
the star. The difficulty is that the gas that does this is also likely to accrete onto the star. Since 
any excess iron sequestered in the doomed planet had to come from the gas in the disk, the gas in 
the disk must be iron poor. If this iron poor gas accretes at roughly the same time as the planet, 
the net change in the surface metallicity of the star will be zero. 

The fact that the net change in metallicity of the star should be small, combined with the poor 
fit of the model to the evolved star data, suggest that accretion must occur after the accretion disk 
vanishes. 

6.2. Ingestion of Gas Giants By Other Means 

Lin (1997) speculates that Jupiter-mass planets might survive the gas disk only to fall onto 
the star later. If this occurs, and if these objects have heavy element abundances similar to gas 
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giants in our solar system, then the star will gain ^ of iron for each accreted object. We can 

find acceptable fits to the data in Fig. (8) if roughly 30% of stars in the solar neighborhood have 
each accreted a single Jupiter-mass bodie after their convective envelopes thinned. These models 
naturally predict that the average metallicity of such stars will decrease by the observed amount 
when their convection zones deepen at the end of their main sequence lifetime. 

However we have more information available to us than just the mean metallicity as a function 
of stellar mass. In particular, we know that the distribution of [Fe/H] is well fit by a gaussion, as 
illustrated in Figure (5). This fact allows us to rule out pollution dominated by late ingestion of 
Jupiter-mass planets having iron contents similar to that of the giant planets in our solar system. 
The idea is to use the sensitivity of the high mass stars to pollution. Adding ~ l.SM^ of iron in 
discrete lumps to 30% of the stars (to give the observed mean [Fe/H]) will produce a population with 
metallicities about three times that of the unpolluted population in these high mass stars. Using 
the Monte Carlo model that reproduces the variation of mean metallicity mentioned in the previous 
paragraph, we obtain the two-peaked histogram shown in Figure (11). This can be compared with 
the observed distribution for those stars above IAMq shown on the same figure. If these stars 
are in fact polluted, the pollution is distributed smoothly over most of the stars rather than being 
concentrated in a fraction of order 30%. We conclude that the metallicity trend we find is not 
produced by the late ingestion of Jupiter-mass objects. 



6.3. Surface Metallicity Enhancement By Mass Loss or Dust Accretion 

The solar wind is observed to be iron rich; in fact many elements with low first-ionization 
potentials are more abundant, while He, with its very high first-ionization potential, is under 
abundant (Meyer 1993; Geiss 1998). If a substantial portion of the stellar convection zone were 
lost in a massive wind then the photospheric iron abundance would be less than the bulk abundance. 
If less massive stars suffer larger mass losses (relative to their more massive convection zones) this 
would produce a population in which less massive stars have lower metallicity than more massive 
stars, as we have found. When these low mass stars evolve, deep convection zones will form and 
mix metal rich material up to the surface; the high mass stars will also show a slight increase in 
iron abundance. 

However, what we appear to see is that as massive stars evolve, their metallicities decrease; 
we don't know what the low mass stars do when they evolve. There is a second difficulty with this 
scenario. The mass loss rates required of low mass stars are excessive; to lose a substantial fraction 
of the convection zone (~ O.OSMq for a O.SMq star) the wind must carry away ~ 5x10^^^ Mq/uv for 
a billion years. The solar mass loss rate is ~ 3 x W~^'^MQ/yr. Recent limits on the mass loss rates in 
K dwarfs are much lower, W"^'^ Mq / yr or less (Lim & White 1996; van den Oord &: Doyle 1997). 

Of course the observation that lithium and beryllium depletions are correlated indicate that 
mass loss does not extend down to masses of order 3 x 1O~^M0 in stars of mass ~ 1.2 — IAMq. 
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Taken together, these facts suggest that mass loss is not responsible for the metallicity variations 
we find. 

Many young stars of roughly solar or larger mass show infrared excesses due to circumstellar 
dust (Aumann et al. 1984; Habing ct al. 1999; Decin et al. 2000). Much of this dust may end 
up on the star. Dust particles which absorb photons moving radially and then radiate them 
isotropically (in the dust rest frame) lose angular momentum in a process called the Poynting- 
Robertson effect. The amount of angular momentum lost {5tj per photon of frequency v is Si/i ~ 
hu/mc^, where m is the rest mass of the dust particle. Thus, for each accreted dust particle, a 
star must emit £/6£ photons, which means it must emit a total energy approximately equivalent to 
the rest mass energy of the accreted material. The total energy emitted during the period of dust 
accretion is larger than this by 1/r, where r is the optical depth to dust absorption. Thus, for a 
star radiating at approximately solar luminosity, the accretion of material will take a time 

l^years f MacA f Lq\ 

— ^\m^)\-L) ^^^^ 

For r ~ 1, this number corresponds roughly to the observed lifetime of dust around young massive 
(F) stars (Habing et al. 1999), but the spherically averaged r is obviously considerably less than 
1; values near 10"*^ are typical. Decin et al. (2000) find longer lived disks (~ 1 Gyr) of similar 
optical depth, but the maximum amount of dust that such stars can accrete is still much less than 
an Earth mass. We conclude that dust accretion is unlikely to provide enough material to explain 
the apparent iron enhancements we see. 



6.4. Accretion of the Interstellar Medium 

Observations of the interstellar medium in the vicinity of the sun indicate that it has roughly 
solar abundances (~ +0.05 dex) of some elements, such as Zn, P, and S (Howk et al. 1999). Iron is 
typically depleted, probably onto dust grains. The youngest stars in our sample have [Fe/H] ~ 0.2 
(Figure 3), slightly larger than the ISM, but our results indicate that about 0.05 — 0.1 of this is due 
to pollution. This is consistent with the notion that the bulk metallicity of young stars will equal 
that of the ISM out of which they form. 

The implication is that older stars have metallicities lower than that of the ISM in which they 
currently reside. Accretion of the relatively metal rich ISM material could in principle produce 
metallicity signatures of the type we have found. However, it seems that the youngest stars in our 
sample may be metal rich compared to the present day ISM. In view of the rather large errors in 
the metallicity determinations we feel that this conclusion has to be considered as tentative. 
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6.5. P Pictoris 

The star /? Pictoris is surrounded by a dusty disk (Smith &: Terrile 1984). The dust is believed 
to have a life time shorter than the estimated age of the star (between 3 x 10^ and 10® years) and so 
must be replenished, possibly by collisions between planetesimals (Lecavelier Des Etangs et al. 1996) 
Transient red-shifted absorption features have been seen in spectra of /3 Pic for the last decade 
(Boggess et al. 1991; Lagrange et al. 1996). These events have been interpreted as infalling comets 
or asteroids. Events with very high velocities (> lOOkm/s) are seen a few times per month 
(Beust et al. 1991; Beust & Morbidelli 2000). Similar transients are seen around numerous other 
young stars (de Winter et al. 1999; Grady et al. 1999). 

The absorption features are believed to be produced by the cometary tails of the infalling 
bodies; Beust Sz Morbidelli argue that the infalling bodies must have radii larger than ~ 10 km 
in order to survive the evaporation long enough to evolve onto the inferred highly eccentric orbits. 
They estimate that l.SM® — ISM® of material is removed from the planetesimal disk by this process 
if it has continued for 10® years. They refer to evaporation of these bodies; however a small but 
substantial fraction of them will survive long enough to strike the star. This could provide enough 
material to give the signature we have seen. 

6.6. Solar Pollution 

The sun is metal rich for stars of its age and mass, as can be seen from Figs. (3) and 
(4). The implication is that the radiative interior of the sun might have a metallicity slightly 
lower than the photosphere. A number of authors (Joss 1974; Christensen-Dalsgaard et al. 1979; 
Levy & Ruzmaikina 1994; Jeffery et al. 1997) suggested this possibility two decades ago as a way 
out of the solar neutrino problem; lowering the metallicity of the solar interior reduces the radia- 
tive gradient, which in turn would lower the inferred temperature in the core. Finally, a lower core 
temperature would predict that fewer neutrinos would be emitted compared to the standard solar 
model. 

The metallicity deficit required to explain the lack of solar neutrinos is much larger than 
the 0.017 dex increase resulting from the accretion of ~ O.4M0 of iron. This avenue for solv- 
ing the neutrino problem was dropped by solar workers, partly for this reason. More recently, 
(Henney &; Ulrich 1998) examined polluted models to see if accurate measurements of the the solar 
five minute oscillations (p-modes) could be used to detect differences between the mctallicities of 
the surface convection zone and the radiative interior. They found that current data was unable to 
distinguish between unpolluted models and models with ~ lOM® of accreted cometary material. 
The effect on the production of solar neutrinos was also negligible. 

As an aside we note that the SOHO satellite has discovered a large number of sun grazing stars, 
a class of comets that plunge into the sun (see (Raymond et al. 1998) and references therein). 
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Roughly 200 have been observed by SOHO over its lifetime to date. Raymond et al. give an 
estimate of about 10^ cm for the radius of the comet they observed. If most of the observed objects 
arc of this size, and if we assume that the currently observed rate has been more or less steady over 
the lifetime of the sun, (one strikes the sun roughly every day as observed in 2000), then the mass 
of accreted material would be about one-one millionth of an Earth mass. Note however that most 
of the comets seen by SOHO appear on dynamical grounds to be the result of the break-up of a 
single object, so that the currently observed flux is likely to be a burst of transient activity. 



7. CONCLUSIONS 

Wc analyzed some 642 stars from the Cayrcl dc Strobel catalogue having spectroscopically 
determined mctallicities, and well determined HIPPARCOS parallaxes. Using a large grid of stellar 
models, we determined the age and mass of the stars in our sample. We then examined the variation 
of metallicity with stellar age and mass. 

We have found striking mass-dependent variations in photospheric iron abundances of main 
sequence solar mass (0.8 — I.8M0) stars in the solar neighborhood. These variations mimic the 
variations seen in lithium abundances. With somewhat less confidence, because of the small sample 
size, we find that the iron abundances of Hertzsprung-gap stars are on average lower than those 
of main sequence stars, and that the metallicities of these slightly evolved stars have no mass 
dependence. These results are consistent with the accretion of an average of OAMg) of iron onto 
the surface of main sequence stars. This strongly suggests that terrestrial- type material is common 
around solar type stars in the solar neighborhood. 
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Planet 




zone boundaries (AU) 


T^Fe {gem 


Mercury 


3.3 X 10-2 


0.22 - 0.56 


1.1 


Venus 


0.29 


0.56 - 0.86 


5.7 


Earth 


0.38 


0.86 - 1.26 


4.0 


Mars 


0.03 


1.26 - 2.0 


0.11 


Asteroids 


0.0005 


2.0 - 3.3 


6.2 X 10" 


Jupiter 


1.8 


3.3 - 7.4 


0.36 


Saturn 


1.7 


7.4 - 14.4 


0.11 


Uranus 


0.9 


14.4 - 24.7 


0.019 


Neptune 


1.0 


24.7 - 35.5 


0.013 



Table 1: Planetary iron content, and the inferred iron surface density of the minimum mass solar 
nebula 



Note. — Data for Mercury, Venus, Earth and the asteroids are from Weidenschilling (1977), for Mars from 
(Longhi et al. 1992), for Jupiter and Saturn from (Guillot 1999b), for Uranus from (Podolak et al. 1991), for Neptune 
from (Hubbard et al. 1995). For the gas giants we take representative values for the heavy element mass; the typical 
uncertainty is roughly ±30%. The zone boundaries are those of Weidenschilling. Using these values, and assuming 
solar abundances, we find the iron masses given in the table. 
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Age (years) 



Fig. 1. — The mass of the convection zone for solar metalhcity ([Fe/H] = 0) stars, plotted as a 
function of stellar age. The stellar mass is used to label the curves. Note that for ages less than 
about 7 million years all stars in the mass range 0.6 — I.GMq have convection zone masses greater 
than 0.1 Mq] stars this young have not yet reached the main sequence. 



-24- 




Fig. 2. — The mass of the surface mixing layer of solar metallicity stars. For masses below about 
I.OMq the mass is assumed to be the mass of the convection zone, as calculated in our stellar models. 
Above 1.5Mq the mass is assumed to be 3 x 10~^Mq, as predicted by meridional circulation models. 
Between these masses the mixed region is chosen by comparison to observations of the lithium dip. 
See the text for a more detailed explanation. 
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Fig. 3. — Stellar mctallicity as a function of the logarithm of the stellar age, where the latter 
is obtained by fitting to our stellar models. Filled squares represent unevolved stars, while open 
triangles represent evolved stars (dwarf stars with surface convection zones ten or more times larger 
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Fig. 4. — Stellar metallicity as a function of stellar mass, where the mass is obtained by fitting to 
our stellar models. As in Figure 3, filled squares represent unevolved stars, while open triangles 
represent evolved stars. 
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Fig. 5. — The histogram of stellar metallicity for 466 main sequence stars (thick line). The best-fit 
gaussian (not shown) has mean ([Fe/H]) = —0.095 and variance a = 0.18. The shaded histogram 
is from a Monte Carlo model of 466 stars. These stars are polluted with a mean mass of O.4M0 of 
iron. See the text for a more detailed description of the model. 



-28- 




Fig. 6. — Average stellar metallicity in age bins of width Alog(a5e) =0.1. There are roughly 50 
stars per bin. The straight line is a least squares fit having slope 0.21 dex/log(Gyr). 
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Fig. 7. — Stellar age plotted against stellar mass. Note that on average the more massive stars are 
younger. Filled squares are actual stars, open squares are Monte Carlo data. It is this relation that 
produces the apparent mass-metallicity trend seen in the unpolluted Monte Carlo model represented 
by the dotted line in Figure 8. 
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Fig. 8. — Average stellar metallicity in mass bins of width O.IMq (filled squares, representing 
unevolved stars). There are roughly 45 such stars per bin. The open triangles represent the evolved 
stars (dwarf stars with surface convection zones ten or more times larger than Mmix- There are 
only 3-6 stars per bin, with a bin width of 0.20Mo for these objects. The evolved stars appear to 
have lower metallicity than unevolved stars of the same mass. The dotted line is the unpolluted 
Monte Carlo model giving the smallest reduced = 3.2. The solid line is the polluted Monte 
Carlo model with an accreted iron mass of O.4M0 and an age-mctallicity slope of —0.14; the fit 
has a reduced = 1.0 The dashed line is the polluted model after the mass of the surface mixing 
region has increased by a factor of ten from its main sequence value, and should be compared to 
the evolved star data represented by the open triangles. This fit also has a reduced of order 
unity. 
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Fig. 9. — The average metallicity, adjusted for age using the slope found in our best fit pohuted 
model (a slope of 0.14). This figure should be compared to plots of lithium abundance showing the 
lithium dip. 
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Fig. 10. — Metallicity versus mass, showing the same data as in Figure 8 and a model in which 
the surface layers of the star are polluted by accretion of enriched material before the star is 10 
million years old. Filled squares are unevolved stars, open triangles evolved stars. The upper 
curve is a polluted model in which a gaussian distributed mass with a mean of lAM^ of iron (c.f. 
the iron content of the giant planets in table 1) is added to the surface mixing zone of a Monte 
Carlo generated sample of stars. For stars lacking convection zones, the mixing zone is taken to be 
3 X 10^'^ Mq. The of the fit is of order unity unity. The lower, dashed, curve gives the metallicity 
of the same sample after the mass of the surface mixing layer has been increased by a factor of ten. 
This gives a reduced = 2.5. 
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Fig. 11. — The histogram of [Fe/H] for stars more massive than 1.4Mq. The thick line is the 
histogram for the observed stars, while the light line histogram, which is also shaded, is the result 
of a Monte Carlo experiment in which a single Jupiter analogue was accreted onto 30% of the 
stars. The bimodal nature of the latter model arises because we have assumed that the planetary 
material, which contains l.SM® of iron is mixed with 3 x 10~^ Mq of the outer layers of the star; 
this results in a large increase in the surface abundance of iron of those stars that accrete the 
planets. The mean metallicity of the entire sample has been forced to match the observed mean 
metallicity. The bimodal distribution of the Monte Carlo model is distinctly different from the 
observed distribution. 



